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ABSTRACT 



The behavior of the B — > niug form factors in the entire physical range 
is examined in a model independent way. Unitarity bounds are further 
constrained by the lattice results (APE). From this analysis, we obtain 
/ + (0) = 0.38 ± .03 if the 5*-pole dominance behavior is assumed. How- 
ever, to get the information on the behavior of the form factor / + , QSR 
results are included. We see the deviations from the pole dominance due 
to the contributions of the higher singularities which are to be treated with 
more precise lattice data. 



1 Introduction 

In this paper we would like to examine the B — > -nlvi form-factors in the entire physical 
region in a model independent way. This is, of course, very important for the forecoming 
B-factories' experiments which will provide us a way to extract the precise value of V u \,. 
This CKM- matrix element plays the crucial role in our understanding of the mechanism 
of CP violation. V u t, is mainly determined from the end-point of the lepton spectrum 
in semileptonic .B-decays. Unfortunately, the theory of the end-point region of the 
lepton spectrum in inclusive B — > X u £i?£ decays is very complicated and suffers from 
large uncertainties. Hence, we explore exclusive decay modes which are also easier for 
experimentalists. The natural candidate is B — ► tt£i?£. However, the determination of V u b 
depends on the knowledge of the physics at large distances, i.e. non-perturbative QCD. 
This problem has proved to be notoriously difficult. We cannot rely on the heavy quark 
symmetry (HQS) to reduce number of the form factors or to fix their normalisationas 
at q^ax- m heavy - heavy transition, heavy quark symmetry helps to calculate the 
normalisation at the zero-recoil point with small and controlled theoretical errors. This 
is the crucial point for the model independent determination of \V c b\. For the case of 
the light resulting particle, we have no such a guidance. The knowledge of the form- 
factor at several points does not mean the knowledge of its functional behavior. Existing 
experimental values (CLEO) [[II] for the branching ratio are model dependent (WBS 
HlOf and Isgur-Wise models): 



The problems with these models are related to the fact that they are always not 
relativistic in some aspects and they do not provide real predictions for the form factors. 
Small recoil behavior of the form factor was also calculated by combined HQS and the 
chiral perturbation theory Ul2| . This approach is model independent, but it is valid only 
in the soft-pion limit. However, there are two methods that are rooted in the QCD from 
the first principles: the lattice QCD (LQCD) and the QCD sum rules (QSR). While 
the lattice can be employed to explore the region close to the zero-recoil point (i.e. 
near g^ax); t ne QSR give us the value of the form- factors for small q 2 . Unfortunately, 
the regions where these two methods apply do not overlap so that the intermediate 
region stays uncovered (at least not with a considerable precision). So far, in the lattice 
approach the values of form-factors were calculated at several points (several q 2 ), and 
then extrapolated according to an ansatz of functional dependence on q 2 . With faster 
computers and new simulations, the extrapolations will be far more constrained. For 
the moment, we would like to obtain the bounds on the form factors, and in a consistent 
way treat the existing results in order to constrain these bounds. In a beautiful series 
of papers, Boyd, Grinstein and Lebed [§, |], H applied an old method || to the heavy 
quark systems. The idea is to use crossing symmetry and dispersion relations in order 
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to relate form factors with QCD perturbative calculations in an unphysical kinematic 
region. To constrain the bounds obtained in this way, we shall use some well estimated 
LQCD-results (APE), as well as some QSR- results. The paper is organised as follows: 
In Sec. 2 we outline the basic theoretical formalism and obtain the bounds. In Sec. 3 we 
strengthen these bounds by using existing results and make an short analysis of the q 2 
behavior of the form-factors. 



2 General Formalism 

For the sake of completeness, in this section we recall the main features of the method 
( see references 0, [5|, 0, ^, ||J]) The current matrix element governing the B — > ir£i>£ 
semileptonic decay is parametrized as 



m.o — mis,. „_l, „m%, — m 2 



< 7r(pO|V(0)|B(p) >= ( P + p'-q " UB 2 "^ r/ + (g 2 )+<f ~ 2 * /V) (2) 



q A q z 



where = vrf^b, and the form factors f + '°{q 2 ) are functions of the momentum transfer 
t = q 2 = (p — p') 2 - They satisfy the kinematical constraint: / + (0) = /°(0). Here, we 
assume that leptons are light so that the physical region of t accessible from this decay 
is < t < (M — m) 2 . The expression for the decay rate is then: 

where A(i) = (t + m 2 B — m\) 2 — 4m|m^ is the usual triangular function. If we want 
to extract the precise value of V U d, it is obvious how important is the knowledge of the 
functional behavior of f + {t). To derive the bounds we consider the two point function: 

W = ij d 4 xe iqx < 0|T(^(x)^ t (0))|0 >= (q»q v - q 2 g^)U T {q 2 ) + g^U L {q 2 ) (4) 

In the case of QCD, the structure functions n T£ (g 2 ) satisfy the once subtracted 
dispersion relations: 

XT,l(Q ) = — 2 l<? 2 =-Q 2 = ~ L U , 02 N2 dt ( 5 ) 



dq 2 nJo (t + Q 

The functions Xt,l(Q 2 ) can be reliably calculated in perturbative QCD as long as we 
stay in the region far from resonances. So, it suffices to calculate them at Q 2 = where 
{m b + m u )AgcD « ( m b + m u) 2 + Q 2 is satisfied. To one loop they read: 



Xt(0) 



37r 2 (m 2 — m^) 5 



2 

{ml - m A u )(m A h + m 4 u - 8m 2 b m 2 u ) - l2m A b m A u log — §■ 



(6) 
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2 

(mf - m 2 u )(ml + m 2 u + m b m u )(ml + m 2 u - Am b m u ) - 6m^m^ log ^ 
Xl(0) = 87i 2 {m 2 h + m 2 u ) 3 " (7) 

For a massless w-quark, it gives: Xt(0) = 1.27 • 10~ 2 /m 2 and Xl(0) = 1.90 • 10~ 4 (from 
now on, we shall simply note xt,l, instead of Xt,l(0)). Also, the 0(a s ) corrections can 
be included [19] and they enhance the above expressions for less than 20%. 



The absorptive parts of the spectral functions ImYlL,T(c[ 2 ) can be obtained by insert- 
ing the on-shell states between the two currents on the l.h.s. of Eq.4. taking the \Bir > 
as the lowest state contributing to the absorptive amplitude. Crossing symmetry states 
that the matrix element is described by the same form factors, but real for t+ < t < oo. 
After integration over the phase-space, the longitudinal part becomes: 

ImU L (q 2 ) > [(t ~^~^ )]1/2 M-l/°(t)l 2 g(t-t + ) (8) 

Below the onset of the i?7r-continuum there is only one resonance - B*(l~), which ac- 
cording to its quantum numbers contributes to It/iIIt, which reads: 

ImU T (q 2 ) > rrf 2 J(t - m|*) + [{t ~ ^ \f + (t)\ 2 e(t - t + ) (9) 

where t± = itlb ± m^. Thus, replacing the absorptive parts in dispersion relations, we 
obtain the following inequalities: 

xl > / Mt)\f°(t)\ 2 dt, 



xt > %+ / Mt)\f + (t)\ 2 dt, (10) 



where we put t~ 2 lmlli(t) = ipi{i)\fi{t)\ 2 6{t—t + ). After integration, the pole-contribution 
becomes very small (f B */m B ,), and can be safely neglected. We absorb this contribution 
intoto xt, anyway. To get informations about the form factors in the physical region, 
we map the complex t-plane onto the unit disc \z\ < 1: 



_ / (m B + m n ) 2 - t 
1 — z V AmBm^ 

By this transformation, the region t_ < t < t + is mapped into the segment of the real 
axis — 1 < z < 0, while the < t < t- is mapped into < z < 1. Two branches of 
the root (or two sides of the cut) are mapped to upper and lower semicircles of \z\ = 1. 
Physically, the relevant kinematic region for the process vacuum — > Bit now lies on the 
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unit circle, while the region for the semileptonic B — > iriug decay lies inside the unit 
circle, on the real axis. Generically rewritten in the z-plane, inequalities (10) are: 

^\M*)fM*<Xi (12) 

2m Jc:\z\=i z 

The functions <f>i(z) are solutions of the Dirichlet's boundary problem @ of finding an 
analytic function on the unit disc. Their values are known on the circle: \<fii{e ie )\ 2 = 
ipi(e l9 ), where tpi are the transformed functions from the integrals (9) . The solution is: 

log 10 (z) | 2 = i I** d6 6 4^ log <p(e* 9 ) (13) 
2ix Jo e — z 

Explicitly, our functions are: 

1 (1 + z) 2 ( ffig + fn-n 1 + z 



MZ 



^GirmBrriK (1 — z) 9 / 2 \2y/m B m w 1 — z J 



v _ m\-m\ 1 + 2: / m B + m w | 1 + z \ 

8mB r m 7T y/2nmBm 7T (1 — z) 5 / 2 \2^/mBm n 1 — z) 

Still, the £>*-pole (z po ze = —0.2519) is below threshold and cannot be ignored when we 
consider / + . We do not know the size of the residue, but the m po i e is known. So, to 
remove the pole, we multiply f + {z) by the Blaschke factor: 

P*(*) = (16) 

1 ^^pole 

which is unimodular on the unit circle and principally does not spoil our analysis al- 
though it will slightly weaken our bounds. Now, the product <f>i(z) f%(z)P*(z) is analytic 
on the unit disc and obeys (12). It should be noted that in our case, there are no branch 
points below the threshold. The last step in deriving the bounds is to construct the 
inner product: 

Let us choose gi(z) = 4>i{z)fi{z) and g%{z) = (1— zz%)~ 1 . From the positivity of the inner 
product, determinant of the (gi,gj) matrix is positive, so that from these two functions 
we have: 

Xi K&WM 

fi(z 2 )(f>i{z2) 



>0, (U 



Wz 2 G IntC 
or explicitly 

X 

[zWl - \z 



m\ 2 < ,^ l21 _Li 2 as) 
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For the case of f°(z) these bounds are depicted on the Fig.l. Coming from general 
principles, the bounds obtained in this way cannot be strong. 




Figure 1: Unitarity bounds on the f°{q 2 ) form factor 
However, we can impose some additional constraints. 



3 Constraints and Analysis 
3.1 Constraints from Lattice data 

To constrain our bounds, first we shall use a set of existing lattice results (for more tech- 
nical details, please see The calculations of f + {t) and f°(t) are performed on the 
APE machine, at f3 — 6.0, on a 18 3 x 64 lattice using the Clover action. The extrapola- 
tion to the bottom mass is done using the scaling behavior predicted by the heavy quark 
effective theory which is fully justified as long as the calculations are performed in the 
vicinity of the zero-recoil point (q 2 ~ Qmax)- The extrapolation to the light quark mass 
is quite smooth and unlikely to be a source of an important uncertainty within a present 
statistical accuracy. For the extrapolations to q 2 = 0, several momenta were injected, 
but the question of this extrapolation left opened though (the results at small q 2 suffer 
from large errors). For our purpose, we take three results with the best precision. 



t = q 2 [GeV 2 } 


17.5 


20.5 


24.4 


26.4 


r 




.69±.10 


.64±.04 


.62±.04 


f + 


1.06±.65 


1.5±.28 


2.47±.26 
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Now, let us define gi(z) = (1 — zz*)^ 1 (i = 1,2,3 lattice points) and again construct 
5x5 matrix whose determinant is positive. Generally, for the case i = 1, . . . , (n — 2), 
this is the n x n matrix: 



1-| Z |2 



l-z„z* 



i—zzi 



1 2 n 2j 



f*(z n )(j)*(z r , 

1 



1 — 22:* 



1-|2„| 2 



> 



(20) 



This is the most compact form of this expression. As it can be noticed on Fig.2, 
the situation for f°(t) is very much improved: while the bare unitarity bounds give 
|/°(0)| < 6.3, constraints imposed by the lattice results concentrated on the opposite 



end (q 2 



Qmax) 



give —0.7 < /°(0) < 2.2 . Errors are included (for the upper(lower) 



bounds we used the upper (lower) limits from the table). 




Figure 2: Constrained unitarity bounds on the f°(q 2 ) form factor; 
3 lattice values are taken (see text) 



The situation for the more interesting f + (t) is not equally good: the lattice errors 
are larger, the Blaschke factor affects slightly ( probably because they apply equally well 
for any value of residue). Still, the bounds are improved. Unfortunately, the analysis 
shows that using the constraint / + (0) = /°(0) does not help for further strengthening 
our bounds. On Fig.3, we show these bounds against the simple pole behavior (dashed 
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line, with / + (0) taken from Ref.fllOl). We see that —0.7 < / + (0) < 2.7 which means 
that the 'bare' allowed range is shrank by a factor of about 4. 




Figure 3: Constrained unitarity bounds on the f + {q 2 ) form factor; 
3 lattice values are taken (see text) 
Dashed line shows WBS model prediction 



Now we can make an analysis of the methods and models applied so far in order to 
see if their predictions of f + (t) behavior fall within bounds in the whole region. For the 
moment, we want to get an information on the function f + (t). 

While LQCD computations give good results around zero-recoil point, the QSR are 
very successefull in the oposite region. To avoid the discussion on error estimations we 
will take some results of the form-factors that are not far from q 2 = point. To relate 
these two sets of results, we will make the fit, but with the essential physics of QCD 
incorporated. Again we perform the conformal mapping: 



1 + z (m B + rriir) 2 - t 



1 - z 



4NmBTn n 



(21) 



which is the same as the previous one for N = 1. Here N is to be adjusted (see below). 
The end points of the kinematic region are mapped as: 



t = 



q 2 

¥max 



rriB + m n — 2y/NmBm 7 , 



Trig + m^ + 2y/Nm B m 7T 

( Vn-i \ 
\Vn + i) 



(22) 
(23) 
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We concentrate on the f + (t). Proceeding the same analysis as above, we obtain: 

(24) 

As long as \z min \ and \z max \ are smaller then one, B — > 7r£z/£ decay posseses a 
small kinematic expansion parameter. Since both f + (z) — > P*(^)/ + (2;) and 0t(z) - ► 
4>t(z) / ^fycr are analytic on the unit disc, we can Taylor expand about z = 0: 

1 oo 

- ^wS-^ (25) 

After the above substitutions, the transformed inequality obtained from dispersion 
relations gives: 



Eqs.(25)and(26) give: 

oo 

E kl 2 < 1 (27) 

n=0 

This is a very important constraint which will be used in what follows. Let us take first 
k-terms of this expansion: 

k is to be chosen in such way that the truncation error is small. The expression for the 
maximum truncation error can be obtained using the Schwarz inequality: 



max \f + (z)- f+(z)\ < max ■ 



E 



y 2n 



yk + l\ 



< max 



P*{z)<t> T {z) 



(29) 
(30) 



We could have applied this method to the previous case (N = 1). However, in this 
case the physical range for B — > -k£u£ ( < q 2 < 2QAGeV 2 ) would have corresponded 
to < z < 0.5183. So, in order to have small truncation error, we should have had 
to include so many terms in the expansion that we could not evaluate such a large 
number of coefficients in expansion with just few points for the fit. However, there is 
a way to disentangle this problem by taking the optimal value of N. In our case, we 
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obtain the optimal value N = 3.45, which gives z min = —0.300 and z max = 0.258. We 
repeat the analysis of the f + {t) constrained by the lattice results and our choice of N. 
To have the small truncation error (±0.03) in the region of available lattice results, 
we have to take k — 4. Beside the three values from the table, we take the fourth 
one / + (19.75) = 1.19 ± 0.24 from the same simulation. The result of this analysis is 
very instructive. In the region around 23 GeV 2 , we see that our bounds are severely 
constrained (Fig.4). 
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Figure 4: More constrained bounds on the f + (q 2 ) form factor; 
3 lattice values and N = 3.45 

It means that we can learn something about the value of the residue for this form 
factor. This point (23 GeV 2 ) is close to the zero-recoil point and is presumably dom- 
inated by the lightest state which couples to the weak current V^, namely the B* 
(M|, = 28.3GeV 2 ). 



(31) 



(F* stands for the residue), we obtain F* = (10.84 ± M)GeV 2 . Of course from 
dispersion relation for this form factor, we have also the contribution of continuum, that 
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we take to be negligible relative to the pole one. If pole- dominance is assumed, this 
value of residue gives / + (0) = 0.38 ± .03. 

3.2 Inclusion of QSR data 

However, our kinematic region is large and we can not that easily conclude about the 
value / + (0),nor the functional dependence of f + {t). To get some additional information, 
we shall take advantage of the knowledge of well estimated values of this form factor in 
the region of the small q 2 . Specifically, we take the updated results of the light-cone sum 
rules which proved to be very stable in the range < q 2 < 15 GeV 2 . In what follows, we 
take only central values in both sets of results ( QSR and LQCD). Again, the goal is to 
learn more about f + (t) behavior. For the fit, we use eight points which ensure that the 
truncation error is less than 0.03. Here, we recall the important constraint (27) which 
eventually reduce the number of coefficients contributing to (28) to four. Completed by 
QSR data, the set for the fit is: 



t = q 2 [GeV 2 } 





4 


8 


10 


15 


17.5 


20.5 


24.4 


z 


0.2584 


0.2239 


0.1829 


0.1590 


0.0856 


0.0380 


-0.0348 


-0.1775 


f + 


0.29 


0.35 


0.45 


0.52 


0.83 


1.06 


1.5 


2.47 



The coefficients a n are: 



a 


<2i 


a-2 


a-3 


CZ4 


3.33-10- 3 


0.36-10" 3 


-33.29- 10" 3 


-5.32-10- 3 


233.11-10 3 



On Fig. 5a, we plot fi(q 2 ) (with constraint (27)) versus f^iq 2 ) (without it). On 
Fig. 5b ft{q 2 ) versus pole behavior, and on Fig.5c f^iq 2 ) versus double pole behavior 
is ploted. For the pole we take (31) with the residue given below the equation, and for 
double pole we take, parameters from Ref . [[Tj| ) 
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Figure 5: a) ft {(f) (fit with data from the table and constraint J2n=o \ a n\ 2 < 1)- Dashed 
curve ftiq 2 ) is corresponding unconstrained fit. b) fi{q 2 ) versus fp i e (q 2 ) (dot-dashed) 
c) ftiq 2 ) versus 'double-pole' behavior (dashed) 



Let us briefly comment the figures. On the first plot, we see the effect of the constraint 
(27). It suppresses the form factor at large q 2 and constrains it to follow rather pole 
dominance. This is evident on the second plot. There is no compelling theoretical 
justification for the usage of the pole- dominance ansatz (except for the points close to the 
zero-recoil). However, from this analysis, we see that this behavior is strongly favored. 
Still, there is no convincing arguments to neglect completely the contributions of the 
singularities above the threshold. According to the present data, we can not estimate 
deviations from the simple pole behavior. Lattice data which are concentrated around 
q^ax should be much more precise for this analysis. There is a permanent tendency 
in reducing all posible sources of errors (statistic and systematic) and we hope that 
next simulations will provide us the data for a pertinent analysis of the excited states 
above the threshold. As we already mentioned, for the extrapolation of the lattice data 
some ansatze are usually considered. These hypothesis must satisfy the kinematical 
constraint (/°(0) = / + (0)), as well as the scaling relations in the heavy quark limit, 
M -> oo W°(q 2 max ~ M ~ 1/2 J + (qmax ~ M 1 / 2 ). Then, the popular ansatz is: 

V) = (32) 

l - -£ 

where n + = n + 1. The UKQCD-resnlts suggested double pole behavior. We take their 
parameters and plot f + (q 2 ) |18| versus 'unconstrained' fit. According to our analysis, 
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the double pole behavior is not acceptable since the constraint (27) bends the curve and 
favours rather the simple-pole one. 

4 Summary 

There is no nonperturbative calculational method to evaluate form-factors in the whole 
kinematic range for the semileptonic heavy to light decays. Even more, there is no fully 
reliable principles which could help. Hence, all that can be done is to use a phenomeno- 
logical ansatz. In this paper we wanted to attack the problem the other way around. 
The calculations performed in the unphysical region are related to the form factors of 
interest via crossing symmetry and dispersion relations. From the derived set of inequal- 
ities, using the conformal mapping, we obtained unitarity bounds on the form factors. 
We have shown that such bounds are not restrictive. To narrow the range of allowed 
values of the form factors, we used results from the simulation on the lattice (APE). 
Additional constraint in this method comes from the optimal choice of parameter N. 
As a result we get a very narrow strip of allowed values of f + (q 2 ) around 23.2 GeV 2 . 
To answer the question about the q 2 behavior of the form-factor, we used the light cone 
QSR results. Having the essential physics incorporated, we perform a simple fit limiting 
ourselves to the central values predicted by the two methods. As a result, we obtained 
the curve which favors behavior dominated by the B*. From / + (23.2 GeV 2 ) we obtain 
the value of residue^] = (10.84±.03)Gel /2 which corresponds, if the pole dominance is 
valid in the whole kinematic range,to / + (0) = 0.38 ± .03. Deviations from this law near 
Imax are present and they come from the excited B* states which are located above the 
threshold. For the proper estimation of their contributions more precise lattice data are 
needed. The other way could be to use the D — > tt£u£ form-factors and to extrapolate 
to B — > ir£ij£ armed by HQS. Then the perturbative contributions must be calculated 
at Q 2 = 16GeV 2 and 0(aj) corrections are large. In our analysis we used the only 
two nonperturbative methods so that the bounds and outlined behavior in q 2 are model 
independent. 
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2 this value of the residue is in a good agreement with the value of qb*b-k obtained in Ref. 
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